Numerical Nahm transform for 2-caloron solutions 
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A new numerical method for performing the Nahm transform for charge k — 2 caloron is presented. 
The Weyl equations with boundary impurities are solved directly and the determination of the 
appropriate basis to the linear system is established. The action densities of the 2-calorons with 10 
moduli parameters are shown. 



I. INTRODUCTION 

Calorons are finite action (anti-)self-dual ((A)SD) so- 
lutions of the Yang-Mills gauge theory on E'^ x S*^ . Along 
with their periodic instanton picture, they can be inter- 
preted as the compound objects of constituent monopoles 
from the perspective of loop group gauge theories P, . 
The intriguing feature of calorons is that they can be 
endowed with a non-trivial holonomy around the peri- 
odic direction, which will serve as an approximation to 
Skyrmions Q . Hence, calorons give a connection between 
instantons, monopoles and Skyrmions. 

There has been formulated a systematic method to 
construct instantons, the Atiyah-Drinfeld-Hitchin-Manin 
(ADHM) construction Nahm has applied the ADHM 
construction for calorons as well as monopoles j5j , known 
as the Nahm construction. In this formulation, gauge 
fields of calorons can be obtained by solving a one- 
dimensional Weyl equation on finite intervals with "im- 
purities" at the boundaries. The dual gauge connection 
in Weyl operators and the impurities are called the bulk 
Nahm data and the boundary Nahm data, respectively. 
The transformation from the Nahm data to the corre- 
sponding gauge connections is called the Nahm trans- 
form. Since the Nahm data are defined in terms of quite 
complicated functions, we need numerical analysis to per- 
form the Nahm transform, in general. 

The Nahm transform for SU(2) caloron of instanton 
charge k = 1 with non-trivial holonomy was, however, 
studied analytically [1, 0|- For the higher charges, some 
exact calorons of fc = 2 were found in @. The authors 
demonstrated that the Nahm data written by the stan- 
dard Jacobi elliptic functions works well for the case of 
calorons. The numerical Nahm transform has been dis- 
cussed in somewhat different context, e.g., (9l.[l0j. 

In this Letter, we perform the numerical Nahm trans- 
form for SU(2) calorons of instanton charge k ~ 2, re- 
ferred to as 2-calorons, in detail. For a given instan- 
ton charge, calorons are classified by the constituent 
monopole charge and mass, respectively. Here we con- 
sider the case of monopole charge (toi,TO2) = (2,2) and 



mass (2/i, /xo — 2fi) = (0,^o), which corresponds to the 
2-caloron without net magnetic charge, and accompanied 
with trivial holonomy. Recently, the 2-caloron Nahm 
data of monopole charge (TOi,m2) — (2,2) of arbitrary 
mass with 16 moduli parameters has been proposed in 
(llj . Since the dimension of framed moduli space of 
SU(2) fc-caloron is 8k [12], this 2-caloron Nahm data 
gives the most general gauge fields. We carry out here 
the Nahm transform for the 10 parameter subset of the 
2-caloron Nahm data whose holonomy is trivial, as a pri- 
mary work for the general cases with non-trivial holon- 
omy. The role of each parameter will be clarified in the 
following sections. 



We should mention the relationship between the earlier 
works on the Nahm transform of calorons and the present 
work. In Q, the authors considered the Nahm transform 
by utilizing the Green function method. In contrast, we 
concentrate here on the Nahm transform by solving nu- 
merically the Weyl equations with impurities. These two 
ways lead, of course, equivalent results for various physi- 
cal quantities such as action density of field configuration. 
However, we expect that the solutions to the Weyl equa- 
tions in the presence of impurities and their numerical 
code considered here make a crucial contribution to the 
progress of Nahm transform for diverse objects jl31] . and 
also D-brane theories, e.g., [l^ . 



This Letter is organized as follows. In section 2, we 
give a brief review on the Nahm transform of calorons. 
In section 3, we make a formulation for numerical anal- 
ysis to the Weyl equations with impurities. In section 4, 
we consider the action density of 2-caloron and give the 
interpretation of moduli parameters. Section 5 is devoted 
for concluding remarks. 
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II. NAHM CONSTRUCTION FOR 
iC-CALORONS 

A. Basic formalism 

In this subsection, we give a brief review of the Nahm 
transform for SU(2) fc-caloron with trivial holonomy. As 
mentioned in Introduction, the Nahm data of calorons 
consists of the bulk Nahm data and the boundary Nahm 
data. The bulk one is k x k Hermite matrices Tfj,{s) pe- 
riodic in s with period fiQ, where /i = 0, 1, 2, 3 and s € I. 
They are smooth functions except for at the boundaries, 
with fundamental interval, say, / — (— /io/2, /io/2). For 
the gauge field to be ASD, they are subject to the Nahm 
equations. 



ds 



T,{s) - i[To(s),T,(s)] - l£,,4T,(s),rfc(s)] = 0, (1) 



where the roman subscripts are 1,2 or 3, together with 
the reality conditions T^{—s) — *T^(s). On the other 
hand, the boundary Nahm data is given by a fc-row vec- 
tor W of quaternion entries, enjoying the matching con- 
ditions 

r,(-Mo/2) - r,(A.o/2) = ^Tr^iW^Wa,), (2) 

where 's are the Pauli matrices and the trace is taken 
over quaternions. 

The link between the Nahm construction of calorons 
and that of monopoles is as follows. In the construction 
of monopoles, the Nahm data should also enjoy ([T]) and 
have simple poles at the boundary of the interval jH, 
[l6| to give correct asymptotic forms of the Higgs fields. 
We can make use of their "non-singular parts" as the 
bulk Nahm data of calorons [l^ , which are piece of the 
monopole data attached with the shorter segment of the 
interval. If we take the monopole limit of the calorons, 
the interval / eventually fills the characteristic period 
of the bulk data so that the simple poles appear at the 
boundaries, accordingly the matching conditions erupt in 
certain cases. 

The caloron gauge fields in the real, or configuration 
space are obtained from the Nahm data through the 
zero modes of the Weyl equations with impurities at the 
boundaries. Denoting a "spinor" U{s; x") defined on / 
with period /xg and a single quaternion V{x°'), the Weyl 
equations are 

{l2fc^ - + x^lk) ® e^}Uis; a;") 



= iW^V(a;")<5(s - Mo/2), (3) 



where 1^ is the kx k unit matrix and = (I2, — icj) are 
the basis of the quaternion. We write these zero modes 
into components as 



where Ui, U2 and Wi, V2 are 2fc-column vectors and 2- 
column vectors, respectively. By using them, we can sep- 
arate dSI into two equations, the bulk Weyl equations on 
/ 

{l2fc^ - + x^lk) ® e^]ui. ^ 0, (5) 

and the matching conditions at the boundary 

Au£ = U£(-//o/2) - (/io/2) = iW'^vi,, (6) 

where i = 1,2. The next step is to find two independent 
pair of the zero modes Wi), (tX2, V2)^ orthonormalized 
as 



j ulubds + vlvb = Sab 



(7) 



where a,b — 1, 2. Putting together these zero modes, we 
obtain the ASD gauge connection of calorons as 



(8) 



As mentioned in Introduction, we have to perform the 
Nahm transform by numerical analysis, which is the main 
aim of this Letter. Note that the numerical Nahm trans- 
form for monopoles was vigorously studied by Houghton 
and Sutcliffe [3. 



B. Exact Nahm data of the 2-caloron 

In this subsection, we introduce an exact Nahm data of 
the 2-caloron, for which we perform the Nahm transform 
in the following sections. In pTj, the bulk Nahm data of 
the 2-caloron on an interval is given in the following form 

Ti{s) ^ fi{s)(Ti -1-51(3)0-3 + dil2, 

12(5)= /2(s)o-2 +^212, 

To = dol2. (9) 

The solutions to the Nahm equations ([T|) are 

/i(s) = a(s) cos 5i(s) = -5(s)sin(/), 

,sn2Ds 

= 

93(3) = ±a(s) sin 0, /3(s) = ±b{s) coscj), (10) 



with 



(a(s),6(s)) = (^D 



k' dn2Ds\ 
- ^ u - 



cn2£'s' cn2Ds / 



, dn2Ds k' 
or D——,D- 



(11) 



cn2Ds ' cn2£'s^ 
where sn, cn, dn are Jacobi elliptic functions of modulus 



f/(s;x") = (mi,M2), V{x°') ^ {vi,V2), (4) /fc, and k' = Vl-k^. The monopole hmit of this bulk 
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data is obtained by putting /xo/2 — >■ 1 and D K{k)/2. 
Note that (f) is not a physical parameter for the trivial 
holonomy calorons and also monopoles, because it can be 
removed by a spatial rotation (see section IV). However, 
we keep it for the purpose of generalization to the non- 
trivial holonomy cases, and also the reliability check of 
the numerical code. 

Next, we consider the matching condition ^ for the 2- 
caloron. In accordance with [ll| , we employ the following 
parameterizations for the boundary Nahm data 

W^=(Al2,pg), g = g^e^, go = cosi/', 

q= (sin -f/; sin sin (p, sin?/) sin cos iy9, sin -0 cos 6*). (12) 

Then ^ reads 

r,(-/io/2) - T,(mo/2) = \pq,(T2- (13) 

From the boundary values of the bulk Nahm data given 
above, we find 9 = 7r/2, = or tt, and 



sn2i?(AW2) . 

± 2Dk — — = Xpsm-th. 

cn2D(^o/2) 



Note that all the parameters are assumed to be real 
number. Consequently, the moduli parameters of the 2- 
caloron of trivial holonomy are the following 10, 

di,d2,d3,do,<p,k,D,X,p,'ip, (15) 

which are subject to the constraint ([T4|. 



III. 



(14) 



NUMERICAL NAHM TRANSFORM: 
FORMULATION 

A. Solving the bulk Weyl equation 



In this section, we give the strategy on finding out the 
numerical solution to the Weyl equations with bound- 
ary impurities, which is based on the construction of 
monopoles a la Houghton and Sutcliffe [l^. From the 
Nahm data (|9]), the bulk Weyl equations ([3]) can be writ- 
ten for the components of ui = ^{uu, U2i, u^e, uu) as 



/ /3 
.93 



1 o\ 

1 ds 

O 1 / V /l - /2 

f X3+ ixo Xi — 1X2 

Xi + 1X2 —X3 + ixo 

X3 + ixQ xi — 1X2 

\ Xi -f 1X2 




ds + ido di — id2 

di + id2 —ds -|- ido 

Q ds + ido di — id2 

di + id2 —ds + ido 



(16) 



This system of ordinary differential equations can be 
solved by the Runge-Kutta method with appropriate ini- 
tial conditions. For the Nahm transform of calorons, 
as well as monopoles, what we need is the basis of 2- 
dimensional vector space spanned by the solutions to 
(fT6|) normalizable on /. To employ the monopole con- 
struction procedure, we have to integrate (|16|) starting 
from the simple poles of the Nahm data. The positions 
of poles are Spoio — ±K{k) /2D, obviously read from 
(fTO|) and (fTT|) . From the Nahm construction of calorons, 
Spoie's are assumed to be located on the outside of /, i.e., 
|spoie| > Aio/2. If we represent pB]) at each pole in the 
following concise form 



(17) 



that the eigenvalues at each pole are, respectively. 



1 

2' 



4±2v^rrT 



(18) 



for an arbitrary </>. Hence, we find the vector space di- 
mensions are three for each pole. In terms of Frobenius 
expansion at the poles, we can find a couple of free pa- 
rameters. Appropriately tuning these parameters, we ob- 
tain the three independent initial conditions at each pole. 
In this way, we can compute the array of independent so- 
lutions u^^\u^^\ ■ ■ ■ ,u^"\ where n is at most 6. 

The 2-dimensional basis u^, {£ = 1, 2), can be obtained 
by taking linear combinations of these solutions. 



Ui = {u 



(1) ,,(2) 



(19) 



where Bs is a regular matrix, then the vector space di- 
mensions of independent solutions which are normaliz- 
able are given by the number of positive eigenvalues of 
the matrix Bg [16|. In the present case, we evaluate 



where u^i are n-column vectors, i.e., u)i 
{lou,uj2i, ■ ■ ■ ,i^ni), to be fixed below. The evaluation 
of U3i is not so straightforward for calorons, as well as 
monopoles. For the monopole construction [l^, uji is 



(ii) Estimate Au'^^\ Au^'^\ Au^'^\ Aw^^) 



(i) Compute from 
both poles 




^Spole -^^0/2 /'-0/2 Spole 



FIG. 1: The method of solving the boundary Weyl equation removing the over-determination. We compute the bulk Weyl 
equation from both poles independently. Then four of them are used to determine linear combinations which can be performed 
with the Gauss elimination for the coefficient matrix of ((24}. 



determined by the condition that a linear combination 
of the solutions, integrated from one pole, matches a lin- 
ear combination of the solutions from the other pole, at 
the center of the interval. In this way, we find the 2- 
dimensional basis of the vector space. For the caloron 
construction, the situation is slightly different. We are 
able to compute the solutions in the whole interval at 
once without taking linear combination. Instead, we have 
to find out the solutions to the boundary Weyl equations 
(|6]) simultaneously, which should be consistent with the 
boundary Nahm data. 



B. Solving the boundary Weyl equation 

For an appropriately defined ue, the boundary Weyl 
equations ^ with the boundary Nahm data (jl2l) are 




A 
A 

p cos ip p sin ip 

-psinip pcosip 



vu 



(20) 



where Ait^ = *(Aui£, Au2f , Am3£, Au4£) and vi = 
^(y\i,vi(). The equations seem to be an over-determined 
system for vii,V2e, but that is not true because the left- 
hand side is not determined at this stage. Our goal is to 
determine the coefficients of the linear combination uji 
in and vi, simultaneously, for the given bulk Nahm 
data • • • , it^"). This procedure can be reduced 

to the problem of linear algebra as follows. 
First, wc solve the upper two rows of (PH)) as 



vn 



-iAM£i/A, ve2 ^ -iAun/X. 



(21) 



By substituting (PT|) into the lower two rows, we find (120)) 
becomes the constraints on Au{ as 

Am£3 — p( Auf 1 cos ip + Aue2 sin t/j^ /A, 

Aum — p(Am£2 cos ip — Auii sin ^/;) /A. (22) 



The next step is to fix a;^ which is consistent with 
From the definition (IT51) . Aug is expanded as 

Aue = (Am(i), Am(2), . . . , Am(")) • uje, 



(23) 



where Am^ = tt^ (-/io/2) - m('Hmo/2). Note that Aug 
and Am^*'''s are 4-column vectors. Now, we can rearrange 
, by using (j23p . into a linear equation for Aui 



p cos ip p sin ip 
—psinip pcosip 



-A 


• (Am(i),Am(2) 





-A 



,Am("))-u;^ =0. (24) 



This can also be regarded as the linear defining equation 
for u^i. Since the number of constraints on a;^ in (j24p is 
two, we have to take at least four solutions of the bulk 
Weyl equation, to obtain two independent solutions 
i.e., the 2-dimensional basis ui. To determine the in- 
dependent solutions u}i to we simply perform the 
Gauss elimination (FiglT]). Having obtained the indepen- 
dent uj^, we take the linear combination (j23l) . and find the 
components of Aug automatically satisfy 1^^. The solu- 
tions of the boundary Weyl equation have already been 
given by ([2T|) . Our method to solve the Weyl equations 
is schematically illustrated in FigH] 



C. Realization of the action density 

To investigate the configuration of calorons in the 
configuration space, it is useful to visualize the action 
density, which is gauge invariant and positive real def- 
inite. With the fact that the field strength of calorons 
Ff_ii, = d^A^ — dyAfj_ -f [Ap, Ai,] is an anti-symmetric ten- 
sor and the (A)SD conditions Fqi — ±^32, -F02 = ±^3i, 
^03 = ±^21 , the action density of calorons can be written 
as 



S{x,t) 



-2Tr2(Fi 



Fk). 



(25) 




FIG. 2: Action isosurface plot of the 2-caloron with /io/2 

Thanks to the (A)SD conditions, we do not have to cal- 
culate the "time" t = xq derivative in the field strength 
so that we can regard t as a parameter in the calculation. 
We perform the visualization of the action density of the 
caloron, mainly by Mathematica (Tsj . 



IV. NUMERICAL NAHM TRANSFORM: 
ANALYSIS 



A. The patchwork 

Combining all these procedures, we now compute the 
action density of the calorons, a typical gauge invariant 
quantity. However, when we carry out the program, we 
always observe unexpected singularities of line form. In 
the calculation of the zero modes of the Weyl equation, 
we always have a trivial phase factor, which depends 
on the configuration space coordinates. These factors 
inevitably cause effect to the gauge field through a fi- 
nite difference method of ([5]). For the analysis of the 
monopoles [l6| , on the contrary, no such singularities oc- 
cur. The finite difference of the zero modes is taken after 
the trace for quaternion, which successfully cancels jump 
or twist of the phase of the zero modes. 

These singularities, however, have no essential ingre- 
dients so we employ the following simple procedure to 
remove them. We can identify numerically the location 
of the singular lines, which depends on the initial condi- 
tions of the Runge-Kutta method. For example, a solu- 
tion has a singular line located on the upper hemisphere 
of the configuration space while another has it on the 
lower hemisphere. Thus, we obtain the regular action 
density by a patchwork with two or three parts of the 
solutions with different parameter sets, by choosing their 
intersection is regular. We employ 50 grid points for the 
dual space, and 61x61x61 lattice points for the con- 
figuration space, which are sufficient for the numerical 
convergence. 
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-- 0.5, t = 0, ^(a;,*) = 2.94, k = 0.0,0.7,0.9 respectively. 

B. The moduli parameters 

Having established the formulation for calculating the 
action density, our main concern is to clarify the role of 
the moduli parameters of the 2-calorons, which are listed 
in (HSl). 

First, we consider a very simple case 

= 1^2 = c?3 — c'o 0, fc = 0, 
D ^ K{k)/2 = TT/4:, p^X, ij^7r/2 (26) 

with the positive sign in ([H]). In this case, A can easily 
be solved as 

where /3 — 2n/fio, which gives the rotation symmetric 
2-caloron of trivial holonomy obtained in [l7| . The con- 
stant surface of the action has a toroidal shape, as ex- 
pected. As the time coordinate t varies, the surface still 
keeps its torus configuration: When we increase t from 
t = 0, the action density gradually reduces its absolute 
value and almost fades away at t = /3/2, and it returns 
to its initial value at t = (3. 

Now, we consider the general case accompanied with 
all the moduli parameters (1151) . The quadruplets di, d2, 
da, do are the parameters of parallel transformation in 
R'^ X 5'-^. Substituting ([9|) into ([3]) one can easily see that 
for a parameter scaling d^i -1-5^, (5^ is absorbed into 

as Xf^^ I y x^ -\- 5^. 

It is already known for the BPS 2-monopoles that the 
modulus fc is a parameter of the geodesic motion of the 
monopoles. The situation is similar for the 2-calorons. 
In Fig. [5J we present the action isosurface plots for sev- 
eral values of A; = 0.0, 0.7, 0.9, which show the dynamical 
motion of the constituent monopoles. 

Next, we consider D, which can be identified as a scale 
parameter of the caloron. Under a scale transformation 
with a being an arbitrary constant, 

Di-^aD, s^s/a, A H> \/aA, p^\/ap, (28) 

the bulk and the boundary Nahm data, (flUl) . ([Ti)) are 
invariant. If we fix the poles of the bulk data located on 





FIG. 3: The action densities on a;iX3-plane with 

/io/2 = 0.5, t ^ 0, ip = Tv/2, for 7 — A/p = FIG. 4: Similar with the Fig[3] but with the case 

3 (top), 7 = 1 (middle), 7 = 1/3 (bottom), re- of ip ^ 37r/20. 

spectively. 



s — ±1, which gives a constraint on Z? as D = K{k)/2, 
then we can find a similar functional form for the bulk 
data as we change the value of D and fio, simultaneously. 
This means that the change of D leads the rescaling of 
the range of/, i.e., ^q. Now, the Weyl equation ([3]) keeps 
invariant if the additional conditions are satisfied 



^ ^fj, 



ax^, U<^U'is/a,x')=y^Uis,x^),{29) 



which obviously describe the spatial rescaling of the so- 



lution. Note that as s reduces, a;^ increases at the same 
order. Hence, we conclude that D works as a scale pa- 
rameter of the solution. 



The triplet. A, p, and ip, are the boundary data of the 
2-caloron. As discussed in [l7| . the symmetric calorons 
are defined such that they are invariant under the group 
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TABLE I: The role of each moduh parameter. 



parameters 



role 



di, d2, ds, do 
k 
D 

X, p, 



Parallel transformation in [ 



X 



Geodesic motion of two BPS monopoles 
Scale of the solution 
Angle of the rotation around a::2-axis 
Parameters of the boundary data 



G C SO (3) as 



(30) 
(31) 



where Or = Rk ® R2 and tr is a unit quaternion. 
R2 G SU(2) is the rotation matrix acts on the quaternion 
coordinate. Similarly, R^ is the /c-dimensional irreducible 
representation of a rotation which acts on fc-column vec- 
tors. One can easily confirm that for rotations about 
X2-axis, (|3ip is satisfied only for X = p,tp = tt/2. This 
means that the calorons are not always axially symmet- 
ric even if the bulk data is invariant. We summarize 
figures of the action density to the non-axially symmet- 
ric 2-calorons with the same bulk data to the rotation 
symmetric 2-calorons in Figs. [3] and SI 

Finally, we consider 0, which describes a rotation 
around a;2-axis in M"^, as mentioned in Section II. This 
is because the bulk data can be written as (for positive 
sign) 





a{s)ai 

/2(s)(T2 
6(s)(T3 



(32) 



and 4> can be set to by the rotation of quaternion a;^e^ 
by i?2, together with the rotation on the Weyl spinors. 
Thus, this parameter is dummy in the case of the 2- 
caloron of trivial holonomy. For the 2-caloron of non- 
trivial holonomy, however, this becomes a crucial param- 



eter because it has a meaning of the relative rotation 
angle of the constituent monopoles. 

Hence, we have found full understandings for all pa- 
rameters. The results are summarized in Table HI 



V. CONCLUSION 

In this Letter, we have performed the Nahm transform 
numerically by directly solving the Weyl equations with 
boundary impurities for the 2-calorons of trivial holon- 
omy, with 10 moduli parameters subject to 1 constraint. 
For the boundary equation, we have proposed the system- 
atic method to extract appropriate 2-dimensional basis of 
the system by taking linear combinations of solutions of 
the Weyl equation. The action density plots are qualita- 
tively similar with those of [Sj], which are based on the 
Green function method. The method presented in this 
Letter is essentially equivalent to [sl]. 

For the most general 2-calorons, with non-trivial holon- 
omy, the program will also work well, though we need to 
solve more complicated matching conditions. The role of 
the moduli parameters in this object is still not fully un- 
derstood. The analysis for such case is now in progress. 

The Atiyah-Manton construction of Skyrmions has al- 
ready been generalized to the case of finite temperature, 
by making use of calorons of trivial holonomy [l^, [l^ . 
For the study of a high energy collision of the nucleon or 
a high density phase of the nuclear matter, it will be valu- 
able to consider the approximation to Skyrmions in more 
general situations. It is extremely interesting whether or 
not the calorons of non-trivial holonomy play a signifi- 
cant role in this subject. The analysis will be reported 
in forthcoming articles. 
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